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Abstract 

We express the Hamiltonian of the quantum trigonometric Calogero-Sutherland model related to the 
Lie algebra D4 in terms of a set of Weyl-invariant variables, namely, the characters of the fundamental 
representations of the Lie algebra. This parametrization allows us to solve for the energy eigenfunctions 
of the theory and to study properties of the system of orthogonal polynomials associated to them such 
, as recurrence relations and generating functions. 

o 

O 1 Introduction 

Q_i< Integrable models play a prominent role in theoretical physics. The reason is not only the direct phe- 
fH ! nomenological interest of some of them, but also the fact that they often provide some deep insights into 
I the mathematical structure of the theories in which they arise. Sometimes, they even reveal unexpected 
relations among different physical or mathematical theories. In classical mechanics, integrability not only 
shows up itself in some of the most important and time-honored problems, such as the Keplerian motion 
or the Lagrange or Kovalevskaya top. It appears also in a plethora of new hypothetical, highly nontrivial 
^ ■ sytems discovered mainly during the three last decades of the past century (see [1, 2] for comprehen- 
sives reviews). Among these, the so-called Calogero-Sutherland models form a distinguised class. The 
first analysis of a system of this kind was performed by Calogero [3] who studied, from the quantum 
standpoint, the dynamics on the infinite line of a set of particles interacting pairwise by rational plus 
quadratic potentials, and found that the problem was exactly solvable. Soon afterwards, Sutherland [4] 
arrived to similar results for the quantum problem on the circle, this time with trigonometric interaction, 
and Moser [5] showed that the classical version of both models enjoyed integrability in the Liouville 
sense. The identification of the general scope of these discoveries comes with the work of Olshanetsky 
and Perelomov [6, 7], who realized that it was possible to associate models of this kind to all the root 
sytems of the simple Lie algebras, and that all these models were integrable, both in the classical and 
in the quantum framework [8, 9]. Nowadays, there is a widespread interest in this type of integrable 
systems, and many mathematical and physical applications for them have been found, see for instance 
[10]. 

The eigenfunctions of the Calogero-Sutherland Hamiltonian associated to the root system of a simple 
Lie algebra L are proportional to some polynomials which form a complete orthogonal system in the 
quantum Hilbert space. For the specials values n a = 1, where g a = K a (n a — 1) are the coupling 
constants, they coincide with the irreducible characters of L. For L = A n , these polynomials provide 
natural generalizations to n variables of the classical orthogonal polynomials in one indeterminate. In 
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particular, for the case with a trigonometric potential, one obtains a generalized system of Gegenbauer 
polynomials. As it was shown in the papers [11, 12, 13], these generalized Gegenbauer polynomials obey 
a set of recurrence relations which constitute a K-deformation of the Clebsch-Gordan series of the algebra. 
The finding of these recurrence relations opened the way to obtain many concrete results on the system of 
polynomials, as for example explicit expressions, ladder operators or generating functions [14, 15]. The 
recurrence relations are also the key ingredient to formulate a perturbative approach to most general 
among the Calogero-Sutherland models, that involving the Weierstrass p-function as potential [16]. 

The aim of this paper is to extend some of the results which have been obtained for A n to the 
polynomials related to other simple algebras. We think that it is a good idea to begin with a concrete 
case. We choose to work in the first place the problem associated to D4 because of the triality symmetry 
exhibited by this algebra, which will help us in simplifying of the treatment. The organization of the 
paper is as follows. In Sect. 2, we explain how to express the Calogero-Sutherland Hamiltonian in terms 
of the fundamental characters of the algebra and how to solve the Schrodinger equation. Then, in Sect. 
3, we obtain the main recurrence relations among the polynomials and use them to give algorithms to 
calculate some subsets of them. Sect. 4 is devoted to find the generating functions for some classes of 
characters and monomials functions of -D4. More recurrrence relations and some other relevant results are 
included in Sect. 5, and finally, in Sect. 6, we give some brief conclusions. Also, we offer two appendices. 
In Appendix A, for the convenience of the reader, we collect some of the basic facts about D4 which we 
use in the main text. In Appendix B we list some polynomials, characters and monomial functions. 



2 The eigenvalue problem 

The Hamiltonian operator for the trigonometric Calogero-Sutherland model related to the root system 
of a simple Lie algebra of rank r has the form 

H = ^(p,p)+ K a (K a -l)sm~ 2 (a,q), (1) 

aeR+ 

where q = (q±, . . . , q r ), p = (pi, . . . ,p r ), ( , ) is the usual euclidean inner product in R r , R + is the set of 
positive roots of the algebra, and n a are constants such that n a = up if || a ||=|| f3 ||. In particular, for 
the case of the algebra D4 (see Appendix A), this leads to the following Schrodinger equation: 

1 / 4 4 \ 4 d 2 

H = --A + «(k-1) MTsin \ qj - qk ) + ^ sin- 2 (q J+ q k ) , A = £ — 5 . (2) 

\j<k j<k J j=i oq J 

The q coordinates are assumed to take values in the [0,7r] interval, and therefore the equation can be 
interpreted as describing the dynamics of a system of four particles moving on the circle. Let us notice 
that there is not translational invariance. We recapitulate some important facts about this model which 
follow from the general structure of the quantum Calogero-Sutherland models related to Lie algebras [9]. 
The ground state energy and the (non-normalized) wavefunction are 

E (k) = 2(p,p)K 2 = 28 k 2 , 

*o(tf) = I II sin fe - Ik) sm( qj + q k ) 1 , (3) 
[j<k J 

where p is the standard Weyl vector, p = \ J2aeR+ a -> with the sum extended over all the positive roots 
of D4. The excited states depend on a four-tuple of quantum numbers m = (7721,7712,771,3,771,4) 
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E m (n) = 2(\ + k P ,\ + k P ), (4) 

where A is the highest weight of the irreducible representation of D4 labelled by m, i.e., A = Ya=i m i^ii 
and Aj are the fundamental weights of -D4. By substitution in (4) of 

*m(?) = *0(9)*m(?), (5) 

we are led to the eigenvalue problem 

-A^S. = £m («)*S. (6) 

with 

A " + "£ ( ctg( * " n) (4 -4) + ctg( * + *> (4 + 4)) (7) 

and 

£m(«) = £?m(«) - ^o(k) = 2(A, A + 2/ep). (8) 
Introducing the inverse Cartan matrix Aj£ = (Aj, A&), we can give a more explicit expression for e m (K): 

4 4 

e m (/c) = 2 ^ A^. 1 mjinfc + 4k ^ A^. 1 = 2(mf + m| + m|) + 4m 2 , + 2 ( m i m 3 + mim 4 + m 3 m 4 ) 

j,k=l j,k=l 

+ 4m2(mi + m3 + 771,4) + 12«(mi + 7713 + 7774) + 20kt772. (9) 

The main problem is to solve equation (6). As it has been shown for the case of the algebra A n 
[11, 12, 13], the best way to do that is to use a set of independent variables which are invariant under 
the Weyl symmetry of the Hamiltonian, namely, the characters of the four fundamental representations 
of the algebra D4. Unfortunately, the expression of these characters in terms of the q- variables (which 
play the role of coordinates on the maximal torus of D4) is not very simple. Denoting the character of 
the irreducible representation of maximal weight Aj as Zj , we find 

4 4 

3=1 3=1 
4 4 4 

Z 2 = X i X 3 + ^i^jT 1 + E X i lx 3 

i<3 i<3 ij 

4 1 1 4 



z 3 = xV h-Vxi, 



Z 4 



i=l Xi X i=l 
_ 1 1 4 

x x r— ! 



where Xj = e 2 ** , and x = yfx\x^x^xl. These expressions make the direct change of variables from qi to 
Zk quite cumbersome. We refrain from trying that approach, and choose an indirect route which has the 
added advantage of being also applicable to other algebras in which the expressions for the characters 
are even more involved. We can infer from (7) the structure of A K when written in the z-variables: 

A K = E a jk (zi) d Zj d Zk + ]T [bf (z t ) + nbf (*)] d Zj . (10) 

j,k=l j=l 
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On the other hand, as it is well known [17], the <£ m are polynomials which, with some precise partial 
ordering for the monomials to be described later, start as follows: 



Therefore, making use of (9), we conclude that 

O'jk(zi) = ZjZk + lower order terms, 



&S r) (*) 



C^Zj + d^ 



0,1. 



(11) 



(12) 



3 v l > 3 ■> 1 3 ' 

Now, to obtain the full expressions for these coefficients, we rely on the fact that, for n = 1, the P m 
polynomial gives the character of the irreducible representation of D 4 with maximal weight J2i=i m i\, 
while for k = the same polynomial is the corresponding symmetric monomial function [9]. Both, 
characters and monomial functions, can be computed by using the information available in the literature 
(see, for instance, the "Reference Chapter" of [18]). In fact, the following short list of polynomials 



P (1) (z) 

-^2,0,0,0 \ z ) 


= z\- 


Z2 ~ 1, 


P {1) (z) 


= Z\Z 2 


- Z3Z4, 


(z) 

VL,0,1,0V*7 


= ZlZ 3 


~ Z 4 , 


(z) 




Z1Z3Z4 + Z 2 


P {0) (z) 

^2,0,0,0^*7 


= *\- 


2z 2 



is all we need to obtain A K . By substituting these polynomials in (6) and using (9), (10), (12) and the 
triality symmetry (that here implies that the final expression for A K should be invariant under permuta- 
tions of the indices 1,3,4), we get enough simple linear algebraic equations to fix all the coefficients. We 
give here only the final result: 

\ A" = (z\ - 2z 2 - s) d 2 Zl + [2 z\ - 4 (z\ + z\ + zfj - 2 z x z 3 z 4 + 8 z 2 ] d\ 2 + (z\ - 2z 2 - s) d\ 

+ {z\ - 2z 2 - 8) d Z4 + (2ziz 2 - 6z 3 z 4 - 8z ± ) d Zl d Z2 + (z 1 z 3 - 8 z 4 ) d Zl d Z3 

+ (ziz A - 8z 3 ) d Zl d Z4 + (2z 2 z 3 - 6z x z 4 - 8z 3 ) d z , 2 d Z3 + (2z 2 z 4 - 6z x z 3 - 8z A ) d Z2 d Z4 

+ (z 3 z A - 8zi) d Z3 8 Z4 + (6k + 1) Zl d zi + [2(5k + 1) z 2 + 8(k - 1)] d Z2 + (6k + 1) z 3 d Z3 

+ (6n + l)z 4 d ZA . (13) 



Once the explicit expression for the operator A K in the z variables is given, the Schrddinger equation 
can be solved iteratively. By direct application of A K to z m = z™ x z™ 2 z™ 3 z™ 4 , we find 

4 



= Sr 



2 a: 



(K)Z 

m E* 

ij€T 



E<^ m ~ ai 

i=l 

m-(2a2+oti+oij) 



jei 



y m-(a 2 +aj) _ 



E 4* 



tj 7 m-(a 2 +ai+aj) 
m 



y m-(ai+2a 2 +a 3 +a4) 



4 E a ™ : 

3 ei 



J m-(ai+2a2+03+04+aj) 
m x 



(14) 



where the sets of indices are / = {1, 3, 4} and T = {13, 14, 34}, and 
a l m = 4m,i(mi - 1), Kn = 12m 2 mj, 



M - 



16 nii rrij, 



dm («) = I6m 2 2 - m 2 - k + E 
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All monomials in A K z m take the form z m M with \x as a positive root. Thus, the polynomial P^ has the 
form 

p&*)= E c ^ m_ ^ (15) 

MGQ+(m) 

where we choose the normalization cq = 1 and, if Q + is the cone of positive roots, 

Q+(m) = {/i G Q+ | z" 1 "' 1 is well defined if z x z 2 z 3 z A = 0} . (16) 

The above-mentioned partial ordering of monomials is given simply by the height of fi, i. e. z" 1 "' 11 > 
Z m-H2 ^ ht(ni) < ht{^2)- From (14), the coefficients c M obey the iterative formula 



e m— /^( K ) ^m( K ) 



(17) 



with 



E a m-(/i-a,) c ft~ai + E ^m-(/i-a 2 -a ;; ) C M-(>2+aj) + E C m-(/i-a 2 -a»-aj) ^-("2+0.+^) 
i=l jel ij&T 

+ 2 a m-(/i-2a2-ai-aj) C M-(2a2+Qi+aj) ^m-(^i-ai-2a2-a3-a4) c /i-(ai+2o2+a3+a4) 

ijGT 

+ ^E a m~(/i-ai-2Q2-a3-04-aj) C ^-( Q i +2a 2 +Q!3+a4+Q!j) • 

Along with the explicit expressions for the roots given in Appendix A, it is suitable for the implementation 
on a symbolic computer program. A list of polynomials obtained through the use of this formula is offered 
in Appendix B. 

3 The structure of the recurrence relations 

As it is well known, all the systems of orthogonal polynomials in one indeterminate z, such that P m (z) = 
z m + • • • satisfy a recursive formula z P m (z) = a m P m+ \(z) + b m P m (z) + c m P m -\{z). In particular, the 
orthogonal polynomials associated to the trigonometric Calogero-Sutherland model for the case of two 
particles and Lie algebra A\ are the classical Gegenbauer polynomials, whose recursive formula is known 
to be 

■W'^-- ^^X'.) ''-- i( ' > 

This formula is reminiscent of the Clebsch-Gordan series for A\. In fact, for re = 1 it reduces exactly 
to this Clebsch-Gordan series: the polynomials are the characters of A\ and the coefficents are equal to 
one. Immediately the question arises about the existence of analogous recurrence relations, i.e., with the 
structure of re-deformations of the corresponding Clebsch-Gordan series, for the polynomials related to 
Calogero-Sutherland models associated to other simple Lie algebras. As it was shown in [11], the answer 
turns out to be in the afirmative for all root systems, but to obtain the expressions for the deformed 
coefficients it is necessary to proceed through a case-by-case analysis. Once the coefficients are known, 
many applications are possible. The aim of this section is to fix the structure of the basic recurrence 
relations for the case of D4 and to give a simple illustration of their use. 

We want to study the formulas for ZjP^(z), i = 1,2,3,4. Therefore, as Pm (z) = zi for rrij = (Sji), 
and the recursive formulas are deformations of the Clebsch-Gordan series, we need to know the weights 
of the irreducible representations whose integral dominant weights are Ai, A2, A3 and A4. For the case of 
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Ai, A3 and A4, these representations have dimension eight. On the other hand, if we act on the highest 
weight with the Weyl group in the way explained in the Appendix A, we obtain eight different weights. 
Thus, these representations include only one orbit of the Weyl group and we are done. For the case of A2, 
the representation has dimension 28 and the orbit of the Weyl group containing A2 has only 24 elements. 
But A2 = «i2, the highest root, and thus this representation is the adjoint one and includes a second 
orbit: the Cartan subalgebra, with four elements of weight zero. Let us summarize. 

• Weights in z x : ±Ai, ±(Ai - A 2 ), ±(A 2 - A 3 - A 4 ), ±(A 3 - A 4 ). 

• Weights in z 2 : ±A 2 , ±(A 2 -2A j ), ±(2A 2 - Ai - A 3 - A 4 ), ±(A 2 + Ai - \j - A fe ), ±(A i + A i - 
A fc ), ±(A 2 - Ai - A 3 - A 4 ), 0, with k £ I. 

• Weights in z 3 : ±A 3 , ±(A 3 - A 2 ), ±(A 2 - Ai - A 4 ), ±(Ai - A 4 ). 

• Weights in z 4 : ±A 4 , ±(A 4 - A 2 ), ±(A 2 - Ai - A 3 ), ±(Ai - A3). 

With these weights, the structure of the recurrence relations results to be as follows: 

^1 ^mi,ro2,™3,m4 (^0 ^mi+l,m2,m3,m4 (^) ^mC^O ^mi- I,m2,m 3 ,m4 (^0 ^'m('^) ^mi+l,m2- I,m3,m4 (^0 

+ c m ( K ) -^mi - l,m 2 +l,m3 ,m 4 ( z ) + ^m( K ) ^mi,m2+l,m3-l,m4-l( z ) 

+ e m( K ) P mi ,m2-l,m3+l,m4+l( Z ) + /m( K ) ^mi,m2,m3+l,m4-l( z ) 

+ 5m( K ) ^mi,m2,m3-l,m 4 +l( Z ) 

z "i Pmi,m2,m :i ,m4 ( z ) = ^mi.mj+l.ma^W + ^hi(k) -Pmi^-l^^t 2 ) + ^(k) ^mi±2,m2Tl,m3,m 4 ( 2 ) 

+ B m (K) 3± Pmx,m2T^,rnz±2,mS. Z ) + ^m(^) 4± Pmx ,m 2 =Fl,TO 3 ,TO 4 ±2 ( z ) 

+ C Y m (/t) ± -F > mi=F i )m2 ±2,m 3 =Fl,m4=Fl( 2: ) + -^m(tt) ± Pm 1 ±l,m 2 ±l,m 3 ^l,m^l i z ) 

+ Dm(K-) 3± -Pmi=Fl,m 2 ±l,m3±l,m4Tl( 2 ) + ^m(^) 4± Pmi^l,m2±l,m 3 ^l,m4±l ( Z ) 

+ E m (n) Pmi^l,m2,m 3 ±l,m4,±l( z ) + ^m(ft) ^mi±l,m2,m 3 Tl,m 4 ±l( z ) 

+ ^ m (fi) 4± -f 5 mi±l,m 2 ,m 3 ±l,m 4 =Fl( Z ) + ^ ? m(^) ± Pm 1 ±l,m2^l,m 3 ±l,m 4 ,±l( z ) 

+ G m (fi) P mi!m2;m3>rn4 {z), 

z 3 Pm\,m2,m 3 ,mA (^0 ^mi,m 2 ,m3+l,m 4 (^) "I - ^m('^) ^"mi,rre 2 ,ra 3 — l,m4 (^0 "I - ^m(^) ^mi ,m 2 — l,m 3 +l,m4 (^) 

+ C m (/t) ^2+1,1713— l,m4 (^) ^m(^) Ani-l,m2+l,m3,tn 4 -l(^) 

+ e m (fv) -P mi +l jm , 2 — l,m 3 ,m4+l( z ) "l~ fm( K ) ^mi+l,m2,m 3 ,m4- l( z ) 

+ 5m( K ) ^rai-l,m 2 ,ra 3 ,ra 4 +l( z ); 

^4 ^mi l m2,m3,m 4 (^) ^mi,m2,ro3,m4+l (^0 "I" ^m(^) ^mi,m2,m3,m4-l(^) ^m( K ) ^mi ,m 2 — l,m 3 ,rri4+l (^0 

+ c m (/t) f mi)m2 +i )m3)m4 _i(-z) + d m {n) f mi _i )m2 +i )m3 _i )m4 (-z) 

+ C m (/t) ^ > mi +l )m2 _l !m3 +l )m4 ('2 ; ) + f m ( K ) -fm.i-l,m 2 ,m 3 +l,m 4 ( Z ) 
+ 5'm( K ) Pmi+l,m,2,m 3 -l,m4, ( z )i 
Where S m (K) 1± Pm 1 ±2,m2^\,m 3 ,mS Z ) means B m(K)' + ^mi+2,m 2 -l,m 3 ,m 4 (*0+-Bm (k) 1 ^mi-2,m2+l,m 3 ,m 4 ( Z )i 

etc, and it is understood that all polynomials involving negative quantum numbers are zero. The re- 
currence relations reflect triality in the fact that not all the coefficients appearing in these formulas are 
independent. There are coincidences upon permutations of the quantum numbers, for instance 

"mi,m 2) m3,m4 ^rri 3 ,m 2 ,mi,m4 ^m4,ra 2 ,m 3 ,mi > (1^) 
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and similarly for 6 m , Cm, dm, 4, /in, 5m- In the same fashion, we have also 

- D mi,m2,rri3,m4 u in3 ) m2,mi ) ra4 -°m4,m2,rri3,mi V ia 7 

and similarly for Dm,Em. 

As an example, let us consider a simple case in which only one of the quantum numbers is nonvanishing, 
namely, 

z l ^m,0,0,o( z ) = ^m+l,0,0,o( z ) + a m{ K ) ^m-1,0,0,0 ( z ) + c m( K ) ^m-l,l,0,o(4 (20) 

where we write a m {n) = 0^ 0,0,0 ( K ) an d c m {n) = 0^0,0,0 Using formulae 

P« (*\ - ~ m m{m-l)[AK 2 + A{m-2)K + (m-l){m-2)] Jm _ 2 m(m-l) 2 
m '°'°'° l j " 1 ~ (m-l + «)(m-l + 3«)(m-2 + K) 1 m-l + K^ 1 Z2 + " ' 

_ m 4 k(k — l)(m - 2 + 2k) m 

m ' 1 '°' otZ] " Zl Z2+ ( m + l + 5K)(m + 2K)(m-l + K) Zl + '"' 

we obtain the coeffficients in (20) 

m(m + 2k) (m — 1 + 4k) (m — 1 + 6k) 



Cm( K ) 



(m — 1 + n){m — 1 + 3k) (m + 3k) (m + 5k) 
m(m — 1 + 2k) 



(m + K){m — 1 + k) 

As a byproduct of triality, we can also write other two recurrence relations with the same coefficients: 

z 3^0,m,o( z ) = P(i,0,m+l,o( z ) + a m( K ) Po,0,m-l,o( z ) + c m( K ) ^l.m-l.oW 

z 4 Po,Ofl,m(. z ) = ^0,0,m+l( z ) + a m( K ) ^0,0,m-l( 2 ) + ^(^ -^l.O,™-!^)' (21) 

The first of these recurrence relations can be used to devise an algorithm for the calculation of the 
polynomials of the form Pm,o,o,o( z ) an d Pm,i,o,o( z )- By multiplying (20) by the differential operator 
A K — £m-i,i,o,o(«), the term involving P^,-i, 1,0,0 cancels. Using the explicit expressions (9), (13), we find 

1 1 + 4k 

^+1,0,0,0 = 4(m + K) [ A *> *ll P m,0,0, (^) " 2(m + K) ^ P m,0,0,o(^) 

m(m + 2k) (m — 1 + 4k) (m — 1 + 6k) 
+ ( m - 1 + K )( m - 1 + 3k) (m + «)(m + 3k) "^AO-oW. 

where, from (13), 

[A",zi] = 4(2;? -2z2-8)d Zl +2{ Zl Z3 -8z 4 )d Z3 +2( Zl z 4 -8z 3 )d Z4 
+ 4 (zjza - 3z 3 ^4 - 4zi) <9 Z2 + 2(6k + 1) z x . 

where, from (13), 

[A K , Zl ] = 4(zf-2z 2 -8)d zl +2{z 1 z 3 -8z 4 )d Z3 +2(z 1 z 4 -8z 3 )d Z4 
+ 4 (z 2 z 2 - 3z 3 z 4 - 4zi) <9 Z2 + 2(6k + 1) Zl . 

Once the polynomials Pm,o,o,o( z ) are known, the recurrence relation (20) provides a formula for each 

-^m,l,0,o( Z ) : 

Cm+l(«0 -Pm,l,0,0 ( z ) = z l ^m+l,0,0,o( z ) ~~ ^m+2,0,0,o( z ) ~~ a m+l(^) ^,,0,0,0 ( z ) • (22) 
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4 Some generating functions 

We present in this section the generating functions for some characters and symmetric monomial func- 
tions. Let us consider first the case of the monomial functions with only one non- vanishing quantum 
number in the form Pmfl,o,o( z )- The generating function for this subset is 



F (t,z)= £* mp i 0) OAo(*)- (23) 



m=0 

In terms of the x variables, the general expression for these monomial functions is 

<Uw = Eh m +^ 7m )' ( 24 ) 

and, in particular, we define Pq^o,o,o( z ) = 8- I n these variables, the computation of Fo(t,x) only requires 
to sum the geometric series: 

*fc*>=g(iir^x)- (25) 

The change to the original z variables can be done by the inspection of the coefficients of the powers of 
t in both the numerator and denominator of this rational expression, with the result 



where 



W,z)=^^, (26) 



N (t, z) = 8 - 7zi t + 6z 2 t 2 - 5(z 3 z A - z{) t 3 + 4(z| + z\ - 2z 2 - 2) i 4 - 3(z 3 z A - zi) t 5 



+ 2z 2 t 6 -z 1 t 7 , 

D(t,z) = l-z 1 t + z 2 t 2 -(z 3 z 4 -z 1 )t 3 + {z 2 + zl-2z2-2)t 4 -(z 3 z 4 -z l )t 5 

+ z 2 t 6 - Zl t 7 + t 8 . (27) 

There is an alternative approach. As the monomial functions are eigenfunctions of A(°) with eigenvalues 
£m,o,o,o(0) = 2m 2 , we have 

-i oo 

A(°)F (M) = E™ 2 rp5, , (4 

m=0 

and, therefore, we can write a differential equation for Fo(t, z): 



^A<°> - (t d t f 



F (t,z) = 0, F (0,z) = 8. (28) 



One can verify by substitution that (26) satisfies this equation. When Fo(t, z) is known, we can easily 
obtain the generating function 

oo 

G (t,z)= £* mp S,o,o(*) (29) 

m=0 

by only recalling (20), which for k = is simply 

z l ^m,0,0,0 ( Z ) = ^m+1,0,0,0 ( Z ) + ^m-1,0,0,0 ( z ) + ^m-l,l,0,o( z )- (^0) 
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This gives 



(31) 



with 

M (t,z) = z 2 -4 + (6zi - 3z 3 z A )t 

+ (-8 - 2 z\ - 10 z 2 - z\ + 4 z\ + 2 ziz 3 z 4 + 4 zf) t 2 

+ (10 zi + 5ziz 2 - 3z\z\ - 4z3Z 4 + Z2Z3Z4 - 3zizj)t 3 

+ (8 z 2 - 4 z 2 + 2 z 2 - Z2Z3 + 4 - z 2 zf ) t 4 

+ (-6 «i - 6 z x z 2 - z z zi + z 2 z 3 z 4 ) t 5 + (8 + 6z 2 + 2 z 2 - zf) t 6 

+ (-10*1 +z lZ2 )t 7 + (A-z 2 )t 8 . 

The computation of the generating functions for the characters an< ^ 1 § oes through similar 
arguments. In this case, the eigenvalues are £ m ,o,o,o(l) = 2m 2 + 12m. Hence, 

00 

F 1 (t, z) = ]T t m PSo,o,o( z ) > W = 1 (32) 

m=0 



is the solution of the equation 

1 F 1 (t,z) = 0, Fi(0,z) = l. (33) 



The Weyl character formula implies that the denominator of F±(t, z) should be the same D(t, z) found 
before. Thus, we try an Ansatz 

and obtain the simple answer 

Ni(t,z) = 1 -t 2 . (35) 

Applying the recurrence relation (20)) we obtain the generating function G±(t,z) for the characters 
P (i) . 

Gi(t, z) = y {z 2 - z 3 z 4 i + (*I + *I - 2z 2 - 1) i 2 - (z 3 z 4 - *i) t 3 + z 2 - 21 1 5 + i 6 } • (36) 

5 More recurrence relations and other results 

In this Section, we give the remaining recurrence relations involving the product of a fundamental char- 
acter times a polynomial with only one non-vanishing quantum number. We also comment the existence 
of some peculiar values for k for which the polynomials associated to some special excited states are 
proportional to integer powers of the fundamental state wavefunction. 

To obtain the mentioned recurrence relations, it is necessary to compute the coefficients of a limited 
number of terms of the polynomials involved. Once the form of these terms is known, we can obtain 
the coefficients in the recurrence relations solving a system of linear algebraic equations. We do not give 
here the full expressions for the coefficients of the required terms, because some of them are too long, 
and only list them: 
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zi z? + AzT~ x za H 


0,771,0,01^/ 




+ £ + (7 z^~ 2 + Z? ziz^" 2 z 3 z 4 + 17 z 1 z^~~ 3 z 3 z 4 




1 


F (zh? 1 ' 2 + z^ 2 zl + z^~ 2 z 2 ) H 


Pi,m,o,o( z ) 




ziz^ + Gzi^J 1-1 ^ + Hz™~ l z 3 z A + • • • , 


Po,m,l,l( z ) 




^Z3^4 + /^1^2 n + ---, 






zP + J^- 2 + ^z^- 2 ^ + ..., 


^m,l,0,o( Z ) 




+ I^ 2 ^ + Nz?~ 1 Z 3 Z4 + M z™ -\ , 


-Pm,0,l,l( 2; ) 




z^4 + AT ^"^2 + z" l+1 + • • • , 


-P"m,l,l( z ) 




^4^4 + Pz^ + Q z x z™- x z z z A + i? (zfzj 1 + + z%z\) + 5 


-^m,0,o( Z ) 




z i z T + T z™ + U z\z™~ 1 z 3 z 4 + W z™ +l H . 



The use of the quantities denoted A to W in the previous formulas in the general structure of the 
recurrence relations give the following results: 



• Formulae of type z 1 P^ mfi (z): 



z l P 0,0,mfi( z ) 
Z l Po,0,0,m( z ) 
z 3 Pmflflfli z ) 
Z 3 -^0A0,m( z ) 
z 4 Pm,0,0fl( z ) 
z 4 Po,0,mfl( z ) 



Plfl,mfi{ z ) + W, 
P lfl,0,m(. Z ) + b r , 
^m,0,l,o( z ) + b r , 
Po,0,l,m( z ) + b r , 

Po,0,m,l( z ) + &m 



(K)-Po,0,m-l,l( z 
( K )-fb,0,l,m-l(^ 
(n)Pm-l,0,0,l{z 
{^)Pl,0,0,m-l( z 
{K)P m -l,0,l,o{z 
(K)Pl,Q, m -l,o{z 



with 

, . m(m — 1 + 4k) 

w) = ( m _i + K )( m + 3K ) • 

• Formulae of type zi P(f m (z): 

z l Po,m,Ofi( z ) = Pl,m,0,o( z ) + dm{ K )Pl,m-l,0,0\z) + e m(«)-P(^ m -l,l,l (-z) 
^-Fftm.O.oW = P 0%,l,o( Z ) +dm(«)-Po,m-l,l,o(^) +em(K)^m-l,0,lW 
^■FolmAoW = P 0,m,0,l( z ) + (i m («)^0,m-l,0,l(^) + em(«)^m-l,l,o(^) 



with 

, . 2m(m + «)(m — 1 + 3k) (m — 1 + 4k) (2m — 1 + 6k) 

m ^ K ' (m - 1 + «)(m - 1 + 2/e)(m + 3K)(2m - 1 + 5K)(2m + 5k) ' 

. . m(m — 1 + 3k) 

6m(K) = (m-l + «)(m + 2K) ' 
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Formulae of type z 2 f£ i0 ,o,o ( z ) : 

*2P£fifl,o(z) = P m,lfiA Z ) + fmW Pm-2,1,0,0(*) + 5m(«) i^_i,o,l,l( z ) + ^(«)^,0,0,oW 
^2 ^0,m,o( Z ) = ^l.m.oW + ^0,l,m-2,o(«) + 5m(«)^ ,m-l,l («) + M«) ^0,m,o( 2 ) 

^2 ^0A0,m(^) = ^l,0,m(^) + /m(«) ^0,1,0,^-2 (*) + 5m(«) ^0,l,m-l(«) + M«) ^O.O.mt^) 



with 



fm( K ) 
0m(«) 
h m (n) 



m(m — l)(m — 2 + 2k) (m + 2k) (m — 1 + 4k) (m — 1 + 5k) 
(m — 2 + /c)(m — 1 + k) 2 (to — 1 + 3k) (m + 3k) (m + 4k) ' 
m(m — 1 + 3k) 



(to — 1 + k)(to + 2k) ' 
4 [-3k 3 + 5k 2 + (6m - 1)k + (to 2 - 1)] 
(to - 1 + k)(1 + 3k) (m + 1 + 5k) 



• Formula for z 2 P^ mfifi {z): 

z 2 Po, m ,0, {z) = Po, m+ l,Ofi{ z ) + km{K) Po,m-l,0,o{ z ) +ft»W ^m-1,1,1 W + fcW ^m-2,1,1 W 
+ r m («) ^m-l,0,o( z ) + ^m-l,2,o( z ) + ^m-l,0,2( z ) + s m(^) ^m,0,o( 2; ) 



with 

k m (K) 

9m(«) 

r m (/«) 



4to(to + k) 2 (to + 2k)(to - 1 + 3k)(to - 1 + 4k) 2 (2to - 1 + 4k)(to - 1 + 5k)(2to - 1 + 6k) 
(to - 1 + k)(to - 1 + 2k) 2 (to + 3k) 2 (to + 4k)(2to - 2 + 5k)(2to - 1 + 5k) 2 (2to + 5k) ' 
m(m - 1 + 2k) 



(to — 1 + k)(to + k) ' 

2to(to - l)(m + k) 2 (to - 2 + 2k)(to - 1 + 3k) 3 (2to - 1 + 6k) 
(to - 2 + k)(to - 1 + k) 2 (to - 1 + 2k) 2 (to + 2k) 2 (2to - 1 + 5k)(2to + 5k) ' 

to(to + k)(to — 1 + 3k) (to — 1 + 4k) 
(to — 1 + k)(to — 1 + 2k) (to + 2k) (to + 3k) ' 

-4t ro («) 



(k + 1)(to - 1 + k)(to + 1 + 4k)(2to - 1 + 5k)(2to + 1 + 5k) 
t m (n) = (-1 + 5m 2 - 4m 4 ) + (2 + 25m - 7m 2 - 40m 3 + 2to 4 )k + (20 - 35m - 123m 2 + 20to 3 )k 2 , 
+ (-22 - 115m + 63to 2 )k 3 + (-19 + 65to)k 4 + 20k 5 . 

Finally, we mention that for k = — ^(n — 1), n € N, the polynomials associated to the dominant 
weight which is n times the Weyl vector p are proportional to a power of the ground state wavefunction, 
namely 

Pj {n ~ 1] = (-l) n 2 12 " J f[ sin(g, - q k ) sinfe + q k ) 

[j<k 

This formula can be verified quite easily by direct application of A~2( n_1 ) in the form (7) to the right- 
hand side: one finds that the Schrodinger equation (6) with the appropriate eigenvalue is satisfied. The 
most convenient way to fix the proportionality constant is by performing an analytic continuation to 
complex qi and considering the region X{ 6 R and x\ 3> x 2 S> £3 S> X4 S> 0. Then, the polynomials are 

-i(n-l) 

dominated by the leading order term, P np 2 ~ z^zl^z^z^, and, on the other hand, using the formulas 
for the fundamental characters displayed in Section 2. one finds z\z 2 zzz^ — xfx 2 X3 and Yl^ <k sm(qj — 
qk) sin(qj + qk) ~ — 2~ 12 :r 3 x 2 ,2;3. This gives the proportionality constant written above. 
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6 Conclusions 



In this paper, we have shown how to solve the Schrodinger equation for the trigonometric Calogero- 
Sutherland model related to the Lie algebra D4 and we have explored some properties of the energy 
eigenf unctions. The main point is that the use of a Weyl-invariant set of variables, the characters of the 
fundamental representations, leads to a formulation of the Schrodinger equation by means of a second 
order differential operator which is simple enough to make feasible a recursive method for the treatment 
of the spectral problem. The eigenfunctions provide a complete system of orthogonal polynomials in four 
variables, and these polynomials obey recurrence relations which are extensions of the Clebsch-Gordan 
series of the algebra. The structure of some of these recurrence relations has been fixed and, for particular 
cases, the coefficients involved have been computed. Also, some generating functions for the polynomials 
with parameter k = 1 and k = have been obtained. These generating functions can give some hints 
about the form of the generating function for general k, see [20]. 
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Appendix A: Summary of results on the Lie algebra D4 

In this appendix, we review some standard facts about the root and weight systems of the Lie algebra 
D4 that the reader could find useful to follow the main text. More extensive and sound treatments of 
these topics can be found in many excellent textbooks, see for instance [18], [19]. 
The most convenient explicit representation of D4 is 




Mj k = Ej )k - E 4+ j t 4 +k , j, k = 1,2,3,4 

Bjk = E jA+k - E kA+j , j, k = 1, 2, 3, 4, j < k 

C jk = E 4+jik - E 4+k:j , j, k = 1,2, 3, 4, j < k 



with (Eij)ki = SikSji. The Cartan subalgebra is 



k i=i ) 
and this confirms that the rank of D 4 is four. The matrix commutators 




[h,Mjk] = (cj-c k )M jk , 
[h,Bj k ] = (cj+c k )Bj k , 
[h,Cj k ] = — (cj + c k )Cj k 



allow us to classify the 24 roots in two groups 

a jk (h) = Cj-c k , j^k, 
a fk( h ) = ±(cj+c fe ), j<k. 
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One can extract the following basis of simple roots 

ol\ = ( 1, -1, 0, ) = a 12 , a 2 



a 3 



= ( 0, 0, 



1, 



-1 



"34, 



= ( 0, 1, -1, ) 
a 4 = ( 0, 0, 1, 1 



"23, 



a 



34' 



where we have given the decomposition of these roots in the basis of H* dual to diag(Mjj), i=l,2,3,4. 
The euclidean relations among the simple roots are 



Thus, the Cartan matrix reads 



(<Xi,ai) 

(a2,«j) 
(ai,atj) 



A = 



0. 



i = 1,2,3,4, 
i = 1,3,4, 
i = 1,3,4. 



/ 2 -1 \ 

-1 2 -1 -1 1 
0-120 

\ -1 2 j 



The positive roots are aij,af-, i < j, and they can be classified by heights as indicated in the table. The 



Height 


Positive roots 


1 


CCi, CK2, CC3, CC4 


2 


"13 = "1 + OL2, «24 = «2 + "3; «24 = "2 + «4 


3 


CK14 = «i + «2 + O3, = «i + «2 + 04, «23 = «2 + CK3 + CK4 


4 


a 13 = a l + a 2 + O3 + «4 


5 


a 12 = a l + 2«2 + «3 + 04 



Table 1: Heights of positive roots. 

Weyl group is easy to describe. The Weyl reflection on the hyperplane in H* orthogonal to the root a 
is s a (v) = v — 2 a. Applying this formula to aij,afj, one readily finds that the most general Weyl 
reflection consists in a permutation of the components of v in the basis plus an even number of changes 
of the signs of these components. This gives \ W\ = 192 for the order of the Weyl group. The fundamental 
weights Afc can be obtained from the equation ccj = X)j=i AjiXj. They are 

Ai = ^(2 ai + 2a 2 + a 3 + a 4 ) = i ( 2, 0, 0, 0), 

A 2 = ^(2a 1 + 4a 2 + 2a3 + 2a 4 ) = - ( 2, 2, 0, 0), 

A 3 = ^(ai + 2Q 2 + 2a 3 + a 4 ) = ^ ( 1, 1, 1, -1), 

A 4 = -(ai + 2a 2 + a 3 + 2a 4 ) = - ( 1, 1, 1, 1), 

and the geometry of the weight system is summarized by the relations 

II Ai || = || A 3 || = || A 4 ||= 1, || A 2 ||= V2, 

(Aj,A 2 ) = l, i = 1,3,4, (\ u \ j ) = - i ij = 1,3,4. 
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The Weyl vector is 

1 4 

p = - ^ = = a = 3ai + 5 ° 2 + 3ct3 + 3 ° 4 = ( 3 ' 2 ' 1 ' )' 

a£R+ j=l 

and the Weyl formula for dimensions applied to the irreducible representation associated to the integral 
dominant weight m = m\\i + m 2 A 2 + + TO4A4 gives 



dimr(m)= TT ( " ; , m ± p) = J- 
Jl + (a, P ) 1440 



with 

4 

P = + 1) JJ(m 2 + wij + 2) (m 2 + + m fe + 3) (mi + m 2 + m 3 + m 4 ) (mi + 2m 2 + m 3 + m 4 ) 

i=l 3 j<k 

where the indices j, k take the values 1, 3, 4. In particular, for the fundamental representations, one finds: 

dimr(Ai) = 8, dimr(A 2 ) = 28, 

dimr(Aa) = 8, dimr(A4) = 8. 

Appendix B: Some polynomials, characters and monomial functions 

We list here all the polynomials, characters and monomial functions with total degree lower or equal to 
three up to triality. 

Polynomials 

A^o,o,o( z ) = z i 

PZoM*) = z *-—n Z2 ~ (l + «)(l + 3«) 

/ n 2 2 2(-l + /e) , 2 2 2n 4(-3 + 5k + 6k 2 + 4k 3 ) 

Z = Z\ Z 3 Z/L - ——, (Zf + Zo+ Za) + -r w 77 rz 2 + 

°' 2 '°' 0V ; 2 1 + k (1 + k)(1 + 2k) v 1 3 47 (1 + k)(1 + 2k)(3 + 5k) 

16(-1 + k)(3 + 10k + 3k 2 + 2k 3 ) 



1 + k)(1 + 2k)(2 + 5k)(3 + 5k) 

3 4(-l + «)(-! + 2/e) 



P^ (z) = ^~YTYn Z ^ + (1 + 2k)(2 + 5k) Zl 

4 

p iAi,oW = z 1 z 3 -——z 4 



1 + 3k 

3 _6 6 12(1 + 2k + 2k 2 ) 

^3,0,0,0 W ~ *i ~ —k XZ2 + (1 + k)(2 + k) ^ 4 " (1 + k)(2 + k)(2 + 3k) Zl 



Po\oA z ) = z l-^- K * wi + (1 + K) 6 (2 + K) &4 + &l + 44) - ( f+t) 2 (t+i) (z?Z2 + Z2zl + Z2 " 4) 

6(10 + 17k + 21k 2 + 10k 3 + 2k 4 ) 2 3(30 + 53k + 4k 2 - 15k 3 + 8k 4 ) 
+ 5(1 + k) 3 (2 + k) Z2 5(1 + k) 4 (2 + k) ZxZzU 
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+ 



12/e(8 + 10k + k 2 + k 3 ) 2 2 2 N 12 ( 30 + 119k + 159 ^ 2 + 124k3 + 80k4 + 24k5 + 4 « 6 ) 

5(1 + k )4 (2 + k) (*1 + Z 3+ Z ±) + 5(1 + k)4(2 + k)(4 + 5k) Z * 

16(-30 + 103k + 440k 2 + 359k 3 + 98k 4 + 86k 5 + 20k 6 + 4k 7 ) 
5(1 + k) 4 (2 + k)(3 + 5k)(4 + 5k) 



Y~.it- / 2 2 2 1 + 3K 4(-l + k)k 2 2 4 / 2 2n 

JW*) = ^ 2 - TT^ 2 " (IT? Z1Z3Z4 + (i + k)2(3 + 5k) Zl + \TTW {Z3 + zJ 

4(9 + 27k + 28k 2 + 16k 3 ) 16(3 + 5k + 2k 3 ) 

-Z2 



(1 + k) 2 (2 + 3k)(3 + 5k) (1 + k) 2 (2 + 3k)(3 + 5k) 

2 i -|- 3^ 2(— 1 + k) 5 

= *i*f ~ ~ (TT^F Z2Z3Z4 " (i + k)(i + 2k) + JTTW (zizi + 

4(-l + k)(9 + 19k + 10k 2 + 4k 3 ) 4(-l + k)(-5 + 2k)(1 + 3k) 

+ (1 + k) 2 (1 + 2k)(4 + 5k) ZlZ2 (1 + k) 2 (1 + 2k)(4 + 5k) * 3 * 4 

8(-9 - 57k - 72k 2 + 28k 3 - 2k 4 + 4k 5 ) 
+ (1 + k) 2 (1 + 2k)(3 + 5k)(4 + 5k) Zl 

nK r . 3 <2 2 x 8(-1 + k) 4(12 + 23k- 11k 2 + 6k 3 ) 

P 1>1>1>0 (*) = z,z 2 z 3 - ^(^4 + z 3 z 4 ) - (i - — - - z 2 z 4 + (1 - 2fg)(2 - - 5fc) ^3 

8(3 - 22k + 4k 2 ) 

-24 



(1 + 2k)(2 + 3k)(3 + 5k)' 
™ , x 4 , 2 2 2x 12 16(1 + 5k) 

P 1>0>1>1 (*) = W 4 - — + ^3 + *4) + (1 + 2 . )(1 + 3k) ^2 + (1 + 2fc)(1 + 3k)S 



Characters 



^1,0,0,0 ( z ) 


= 21 




^0,1,0,0 ( z ) 


= z 2 




^2,0,0,0 ( Z ) 


= z\- z 2 -\ 




^0,2^,0,0 (*0 


= z\ + Z 2 - Z1Z3Z4 






= ZiZ 2 - Z 3 Z 4 




P (1) (z) 


= Z1Z3 - 2:4 




P (1) (z) 


= z 3 — 2ziz 2 + 2:32:4 - 2zi 




^0,3,0,0 (^) 


= z 2 + 3z 2 + 3z 2 - 2212:22:32:4 + z 2 zf + z 2 zf + zf z 2 




- (z\ + zf + zf )z 2 - Z1Z3Z4 - Z 2 - z| - 


-z 2 + l 


"^2,1,0,0 ( z ) 


= z 2 z 2 - z| - zi 2:32:4 + zf + z\ — 2z 2 - 


1 


^1,2,0,0 (^) 


= Z\z\ - Z4Z3Z4 - Z 2 Z 3 Z 4 + Zl(zf + zf) 


- Zl 




= ZlZ 2 Z 3 + Z\Z 3 - (zf + Z 3 )z 4 + Z 4 




<1» 


= Z1Z3Z4 - z\ - z 3 - zf + z 2 + 2 





Monomial functions 

,(o) 



^1, 0,0,0 ( z ) — z i 
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4i,o,o(*) = ^2-4 

P 2%o,o( z ) = *i-2*2 

p o%o,o( z ) = z 2 - 2 21*3*4 + 2 + 2 z| + 2 - 4 z 2 - 8 

^M^oW = Z !Z 2 - 3z 3 z 4 + 2zi 

p l%i,o( z ) = z 1 z 3 -Az A 

P 3%o,o( z ) = z i - 3z i z 2 + 3z 3 z 4 - Szi 

P o%o,o( z ) = z% + 6z% + 9z 2 - 3Z1Z2Z3Z4 + 3zlz% + 3zfzl + 3z%zl 

- 3(zj + z\ + zf)z2 - 9Z1Z3Z4 - 4 

P 2%,o( Z ) = Z l Z 2 ~ 14 ~ Z 1 Z 3 Z 4 + 4z§ + 4zf - 6z 2 - 8 

A^O.oC- 2 ) = Z 1 Z 2 ~ 2zfz 3 Z4 - Z2Z3Z4 + 2zf + 5zi(z% + z\) - 9Z\Z2 - 5Z3Z4 - 6zi 

P i°llfl( Z ) = Z1Z2Z3 + 8Z1Z3 - 3(zf + Z^)Z4 + 4z 2 Z4 - Az^ 

p i%i,i ( z ) = z ^ z ± - Az l - 44 - 4^4 + 12z 2 + 16 
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